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O ■ ABSTRACT 
(N 



Using the Gemini infrared camera on the 3-meter Shane telescope at Lick 
Observatory, we have searched for broad-band J and K' photometric variabil- 
ity for a sample of 15 L and T-type brown dwarfs, including 7 suspected spec- 
tral binaries. Four of the dwarfs— 2MASS J0939-2448, 2MASS J1416-I-1348A, 
^ '. 2MASS J1711+2232, and 2MASS J2139+0220— exhibit statistically significant 

. variations over timescales ranging from ~0.5 hr to 6 days. Our detection of 

Oh. variability in 2MASS J2139-I-0220 confirms that reported by Radigan et al., and 

O , periodogram and phase dispersion minimization analysis also confirms a variabil- 

. ity period of approximately 7.6±0.2 hours. Remarkably, two of the four variables 

^ are known or candidate binary systems, including 2MASS J2139-I-0220, for which 

we find only marginal evidence of radial velocity variation over the course of a 
>• . year. This result suggests that some spectral binary candidates may appear as 

such due to the blending of cloudy and non-cloudy regions in a single "patchy" 
lO , atmosphere. Our results are consistent with an overall variability fraction of 

35 ± 5%, with no clear evidence of greater variability among brown dwarfs at the 
O ■ L dwarf/T dwarf transition. 

cn 



1. Introduction 

Low-temperature, low-mass stars and brown dwarfs of spectral types L and T exhibit 
a broad diversity in color and spectral morphology arising from variations in effective tem- 
perature (Te//), surface gravity, composition, atmospheric dynamics and cloud properties 
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( iKirkpatrickl |2005| and references therein). Clouds are a particularly unique and important 
component of L and (some) T dwarf atmospheres, where the photospheric cond i tions are con- 



ducive to the forma tion of mineral and metal condensates (iTsuji et al.lll996l : iHelling et al. 
2OOII : iLodderd |2002| ). Evidence for condensate formation is found in pronounced changes 
in spectral features (e .g., the depletion of TiO and VO through t he L dwarf sequence 
Kirkpatrick et al.l 120001). red near-infr ared colors {J — K ^ 1.5-2.5; IChabrier et al.l 12000 



Leggett et al. 



2001 



Knapp et al.ll2004l) and the direct detection of silicate grain absorption a t 



mid- infrared wavelengths (ICushing et al.ll2008t iBurgasser et al.ll2008al : iLooper et al. 



2008b|). 



While one- dimensional atmospheric models have had modest s uccess in reproducing 



some 



of the spectral and color chara c teristics of cloudy L dwarfs (e.g.. lAckerman fc Marleyll2001 
Cooper et al.l l2003l : iTsujil l2005l : iHelling et al.l 120081 ) , they provide no insight into the distri- 
bution of cloud structures across the surfaces of a rotating brown dwarfs (i.e., the presence of 
banding, vortices and jets), nor can they predict the long-term evolution of these structures. 
Cloud effects are also important for young and/or closely-orbiting hot exoplanets, which have 
comparable equilibriurn temperatures and evidence of unusually thi c k photospheric clou ds 
(iBurgasser et al.l l2010d : ICurrie et al.l 1201 ll : iMadhusudhan et al.l I2OIII : iBarman et al.l I2OIII ) 



While the relatively small sizes {R ~ RjupUer) and large distances of brown dwarfs 
prohibits resolved imaging of cloud featur es, photometric vari ability proves to be an effec- 
tive indirect probe (see recent reviews by iGoldmanI |2005| and lBailer-Jonesll2005l ). Studies 
over the last decade have sh own that perhaps as ma ny as 70% of L dwarfs are variable in 



broad-band photon aetry (e.g., iRockenfeller et al.lbooel ). man y of which are too cool to exhibit 



magnetic spotting ( IGelino et al.ll2002l : iMohanty et al 



2OO2I ). The variatio ns are typically pe- 



riodic over timescales of ~2-10 hours, consistent with rotation periods (IBailer- Jones! 12004 



Reiners fc Basril 12008). There is also evidence of changes in light curv e behavior over mul 



tiple rotation cycles (IBailer- Jones fc Mundtll2001t lArtigau et al.ll2009l ). Variability appears 



to be more common for dwarfs at the transition between the L and T dwarf classes (L8 to 
T5), sources which also exhibit a remarkable brightening between 1-1.3/im ( iTinney et al. 



Bureasser et al. 


2006b: 


Liu et al. 


2006; 


Looper et al. 


2008a) 



sphere ( 


Knapp et al. 


Bureasser et al. 


2002) 



20041 ) or fragmenting of the cloud layer ( lAckerman fc Marleyl 12001 



Unfortunately, the majority of brown dwarf variable detections h ave small amphtudes 



(<0.1 mag), and fe w sources exhibit truly significant light curves ( lArtigau et al.l 12009 



Radigan et al.l |2012| ). This makes it difficult to explore cloud meteorology in a statisti- 
cally robust manner, or quantify the dependence of variability on other parameters such 
as rotation, temperature, surface gravity, metallicity, or overall cloud content. In an effort 
to increase the sample of variable dwarfs and more fully explore brown dwarf cloud mor- 
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phologies, we have conducted a near-infrared variability study of 15 L and T-type dwarfs 
using the Gemini near-infrared camera at Lick Observatory. In Section 2 we describe our 
observation methods, data reduction procedures and photometric analysis. In Section 3 we 
present the light curves for all of our sources, and describe the statistical analyses used to 
identify significant variability and search for periodicities. We discuss individual variable 
sources in detail in Section 4, where we also present additiona l spectroscopic obser vations of 
the previously-identified variable dwarf 2MASS J2139+O22C0 faadigan et"aDl2012h . Results 
are discussed in Section 5 and summarized in Section 6. 



2. Observations 



2.1. Sample Selection 

We selected sources for our variability study using the online DwarfArchives cataloj^ of 
over 800 known L and T dwarfs. We narrowed this sample to sources visible in the northern 
hemisphere (5 > —20°), and focused on two categories of targets. First, we selected L and 
T dwarfs with J < 15.5 and < 15, whose 2MASS J — Kg col ors are >0.25 mag re dder or 
bluer than the mean color for their spectral class as quantified in lFahertv et al. N2OO9I). Such 



sources may have with unusually thick or th in clouds, respectively (e.g., iKnapp et al.l 12004 : 
Burgasser et al.ll2008al : iLooper et al.ll2008bl ). allowing us to examine variabihty dependence 
on total cloud content. We excluded sources spectroscopically known to be young /low-gravity 



objects, which are red d ue to reduced collision-induced H2 opacity (e.g, iKirkpatrick et al. 



2006 



Al 



ers et al.l 12007); and metal-poor subdwarf s, which are blue due to enhanced H2 



opacity ( ISaumon et al.l Il994j : iBurgasser et al.l 120031 ). We required at least two additional. 



comparably bright stars be present within Gemini's field of view (FOV). Our second sample 
comprised spectral bin ary candidates, sources whose spectra appear to be a blend of two 
distinct spectral types (lBurgasserll2007l ). with the goal of measuring potential echpse events 
by these unresolved systems. For these sources we relaxed the declination, brightness and 
color constraints, but not the requirement for at least two bright calibrator stars. The final 
list of targets, which includes four red dwarfs, two blue dwarfs and nine spectral binary 
candidates, is provided in Table [TJ 



^Throughout the text, we shorten source names to 2MASS/SDSS JkhmmzLddmrn, were hhmm and 
±ddmm are the Right Ascension (hours and minutes) and dechnation (degrees and arcminutes), respectively, 
at equinox J2000. Full source coordinates are given in Table [TJ 

^http : //DwarfArchives . or^ 
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2.2. Lick/Gemini Near-Infrared Imaging 



Observations were undertaken at t 



le Lick Observatory 3m Shane Telescope witli tlie 



Gemini dual-channel infrared camera (IMcLean et al.l Il993[ ) . The Gemini camera allows 



simultaneous imaging in two near-infrared bands over a 2'.9 x 2'9 FOV with resolution 
0'.'68 pixel"^, permitting measurement of both brightness and color variability. We used 
the broad-band J (Ac = 1.25 /im) and K' (Ac = 2.12 /im) bands. The former samples a re- 
gion with low molecular gas opacity, and is thus more sensitive to variations in cloud opacity 



( lAckerman fc Marleyll200ll ): the latter samples a region dominated by collision- induced H2 



absorption. 

Data were obtained during two runs, 2010 March 23-29 and 2010 October 25 -November 
1 (UT); each set of observations will hereby be referred to using their respective months (i.e. 
March and October). Observations are summarized in Table |2j For the March observations, 
we performed continuous monitoring over 1-2 hours, obtaining 45s exposures with a repeated 
4-position dithering sequence. For the October observations, we obtained single sets of 5- 
position, 30s dithered exposures interleaved over 2 to 4 hour periods. For each revisit, care 
was taken in repositioning the source on the same pixel at the start of each dither sequence 
to minimize pixel sensitivity variations. 

Data were reduced using custom routines developed in the IDlfl environment. Raw 
image frames were first corrected for instrumental and telescope response by dividing through 
by a normalized fiat-field frame, constructed from median-combined observations of the 
twilight sky typically taken on the same night or one night prior/afterwards. Flat-fielded 
images were dark-subtracted and bad pixels replaced with nearest neighbors. The dithered 
science images were then pair-wise subtracted, manually taking into account source overlap 
in the selection of subtraction pairs. Analysis (Section 3.1) was based on these reduced 
images. 



2.3. Keck/NIRSPEC Near-Infrared Spectroscopy 



One of our targets, J2 139+0220, was identified by lBurgasser et al.l (l2010al ) as a spectral 
binary candidate based on its unusually st rong CH4 abso r ption at 1.2 /xm, and inferred 
component types of L8.5±0.7 and T4.5±1.5. iRadigan et al.l ( 12012| ) report that the source is 
unresolved in Hubble Space Telescope Near-Infrared Camera and Multi-Object Spectrometer 
(NICMOS: lThompsonlll992l ) observations, ruling out the presence of binary companion wider 



•^Interactive Data Language 
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than ~1 AU. To examine whether J2139+0220 is simply a tighter binary, we obtained three 
epochs of high resolution near-infrared sp ectroscopy of this so urce using the NIRSPEC echelle 



spectrograph on the Keck II telescope (IMcLean et al.lll998l ). Data were obtained on 2010 



November 26, 2011 July 6 and 2011 September 7 (UT); in all cases conditions were clear 
to light cirrus, with seeing of 0'/7-l" at J-band. We used the N7 filter and 0'.'432xl2" slit 
to obtain 2.00-2.39 /im spectra over orders 32-38 with A/AA = 20,000 (Av = 15 km s"^) 
and dispersion of 0.315 A pixel"^. An AB nodding sequence was obtained with individual 
exposure times of 1000 s each. We also observed the nearby AO V star HD 198070 (V = 6.37) 
for flux calibration and telluric correction, and dark, quartz lamp and NeArXeKr arc lamp 
frames at the beginning or end of each night without moving the instrument configuration. 

Data were reduced using the REDSPEC packag^. We focused our analysis on order 3 3 



(2.291-2.326 /xm), which samples the strong CO band around 2.3 /im ( iBlake et al.ll2008bl ). 
Image rectification and wavelength calibration (second order polynomial) were performed 
using the arc lamp images, and for each observation, wavelength solutions were corrected 
for barycentric motion. Flux calibration and telluric absorption corrections were determined 
from the AO V spectra assuming a 9480 K blackbody. Figure [T] displays the three epochs of 
NIRSPEC data for J2139+0220, corrected for barycentric motion. Signal-to-noise (S/N) for 
these data are roughly 30-40, and show strong molecular absorption arising from ^^CO 2-0 
and 3-1 transitions. Radial velocities based on these observations are discussed in Section 4.5. 



3. Photometric Monitoring Observations 
3.1. Aperture Photometry 



Aperture photometry was measured for all visible sources in each of the data frames. 
Centroids were found using the ce ntroid routine, a variant of the find routine for the 
DAOPHOT package ( IStetsonI 119871 ). Measurements were obtained in two passes. First, 
sources were located in the first image of a given sequence, and their pixel positions were 
used as initial guesses for subsequent images based on the dithering sequence, allowing for 
a drift of up to 5 pixels. In some cases, offsets between images were adjusted manually. 
These centroids were used to shift, stack and median-combine the full sequence of observa- 
tions to obtain a high signal-to- noise (S/N) master image. An example master image for 
2MASS J2139+0220 is shown in Figure [21 with the source and comparison star marked. In a 
second pass, the master image was used to identify dimmer sources, and positions for these 



http : / / www2 . keck . hawaii . edu/ inst/nirspec/redspec/ index . html 
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were determined by their relative offset from the brighter stars. Aperture photometry was 
then performed, samphng aperture radii of 2-14 pixels (1.3" - 9.4") and a median background 
sampled in 10-30 pixel (6.7" - 20.1") annuli. We determined an optimal aperture for a given 
observation sequence (applied to all measurements of all sources in the sequence) as that 
which minimized the measurement uncertainty: 

4- = + ^ [(^''-^)' + (^^)'] ' (1) 

where F^nt is the integrated counts (DN) in the aperture, corresponding to shot noise; A is 
the fractional number of pixels in the aperture, aback is the standard deviation in background 
counts, and RN = 2.1 DN is the read noise of the detector. We ultimately adopted a 4-pixel 
aperture, which was found to minimize uncertainties for nearly all of the observations in 
both datasets. We then averaged the count measurements (equivalent to integrated fluxes) 
within a dither sequence using uncertainty-weighting. We also computed the average seeing 
in each image from the widths of gaussian fits to each source point spread function (PSF). 



3.2. Differential Photometry and Variability 

For each detected source, we first computed normalized fluxes, defined here as the 
individual flux measurement i divided by the time-averaged mean over a given observing 
sequence: 

F, = F,/{F), (2) 

where (F) = J2iLi ^i- Differential magnitudes were computed for each science target 
(UCD) by dividing its normalized flux by the sum of normalized fluxes of all comparison 
stars (C) in the Gemini FOV: 



UCD 



A = -2.5/05(10 , ^ 



Fr, 



(3) 



Uncertainties were computed as 



4 



1.179 



+ 



N /- n2 



UCD , 



(4) 



where a = a/{F). iBlake et al.l (j2008al ) utilized a similar equation for computing their 
differential magnitudes and uncertainties, but without normalized fluxes. The comparison 
stars were initially chosen as all detected sources in the FOV excluding the science target. 
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Typical uncertainties in differential magnitudes ranged from 1% to of 5% for most 
sources. We conservatively introduce a 1% error in quadrature with these values to account 
for systematic uncertainties associated with, e.g., flat fielding, OH sky line variations, and 
other second-order atmospheric opacity effects. 



quence: 



We then measured the reduced of each source's light curve over the observing se- 



1=1 



To distinguish truly variable sources, both among the science targets and comparison stars, 
we determined a threshold (Xr)mm to exclude false positives, defined here as where the 
highest 1% of Xr values for a non- variable light curve with noise a exceeded the lowest 
10% (or 1%) of values for a variable light curve with semi-amplitude 5*0". Distributions 
of Xr values were generated through simulated data via Monte Carlo methods, assuming 
gaussian noise, 1 < S* < 10, and sampling light curves identically as our observations. We 
found that a minimum threshold was reached only for 5 > 3, which defines our minimum 
detection threshold for each target in units of average uncertainty, and that ixDmin > 1-5 
(2.0) corresponds to a bona-fide variable source at the 90% (99%) confidence level. The 
Monte-Carlo simulation to demonstrate the thresholds is shown in Figure |H The 90% level 
for the non-variable signal is Xr = 1-5 and the overlapping point with the bottom 1% of the 
variable and the top 1% of the non-variable signals is approximately Xr = 2.27. 

In addition to comparing xf- to {xDmin for all science targets, we checked that all 
comparison stars had x^ < 1.5; any exceeding this threshold were rejected (to a minimum 
of 2 sources), and differential magnitudes recomputed. 

We also checked for trends in differential magnitudes with airmass or seeing variations. 
As the spectral energy distributions of our sources are generally distinct from background 
stars, wavelength-dependent changes in telluric opacit y could potentially give rise t o false de- 



tections of variability in broad-band photometry (e.g.. iBailer- Jones fc Lammll2003l ) . Figure [3] 
plots differential magnitudes (minus sequence mean) versus estimated seeing and airmass in 
J- and K-band. In the J-band we found correlation coefficients of 0.007±0.009 versus seeing 
and 0.003±0.004 versus airmass for all dwarfs, indicating that any variation in telluric opac- 
ity was spectrally grey. For the L-dwarfs, we found correlation coefficients of 0.004±0.005 
(seeing) and 0.008±0.010 (airmass), and for the T-dwarfs we found coefficients of 0.02±0.03 
(seeing) and 0.005±0.006 (airmass), again indicating no statistically significant trends. Sim- 
ilar results were found for the K' data. 
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3.3. Results: Variable Sources 

Table [3] and table H] list the reduced Xr values, median photometric uncertainties and 
variability amplitudes of our targets for J and i^-band data. Figures [5] through |6] display 
the differential light curves. Based on the criteria established in Section 2.4, adopting a 90% 
detection threshold {Xr > l-^)' imposing Xr of the comparison sources to be less than 
1.5, we find that J1315-2649 , J1711+2232, and J2139+0220 show significant variability at 
J during one or more observing sequences, while J0939-2448 exhibits variability at K'. We 
also analyzed combined light curves for sources observed over more than one day. In this 
case, J2139+0220 was found to be variable over a period of 3 days. Among our variable 
sources, J2139+0220 and J0939-2448 are the only sources with S > 3, with J1711+2232 a 
marginal case with 5* > 2.5. Each of these sources is examined individually in Section 4. 



3.4. Lomb-Scargle Periodogram Analysis 

For variable sources monitored over periods greater 1 hou r, we perform ed Lomb-Scargle 



periodogram analyses to search for significant periodicities ( IScargld Il982l ). We examined 
frequencies corresponding to periods ranging from three times the sample rate (i.e., above 
Nyquist sampling) to one-half the total observation period. Frequencies were oversampled by 
a factor of ten. To properly account for uncertainties, we calculated periodograms for 1000 
realizations of the data with each differential magnitude adjusted by a normal distribution 
scaled to the corresponding uncertainty. The final "raw" periodogram and its uncertainty 
was taken as the mean and standard deviation in each frequency bin. 

Given the large time gaps in data taken over several nights (e.g., J2 139+0220), sampling 
effects were a concern. To separate false periods arising from sampling from true periodic 
variations, a window function was computed by randomly shuffling the observed light curve 
1000 times and computing the associated periodograms. For each frequency bin we deter- 
mined the maximum value encompassing to 95% of the amplitude distribution. Dividing the 
raw periodogram by this window function yields a significance spectrum. Based on experi- 
ments with simulated light curves sampled in the same manner as the data, we determined 
that peaks in the significance spectrum > 1 correspond to real variations. These peaks were 
fit to Gaussian function, whose peak and half-width were used estimates of the variability 
frequency and its uncertainty. 
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3.5. Phase Dispersion Minimization Analysis 

In addition to periodogram analysis, we performed p hase dispersi o n min imization (PDM) 



analyses on our variable sources using the formalism of IStellingwerj (jl978l ). This technique 
compares the overall variance (a^) of the light curve to the piece-wise variance of the same 
curve folded over a given period (s^). Significant periods appear as a minimum in the ratio 
= s^/(T^. We sampled periods over the same range as the periodogram analysis in steps 
of 0.01 hours, and light curves were sampled in steps of 0.5 in phase. We incorporated un- 
certainties in the measurements in a similar manner as the periodogram analysis, and also 
examined sampling and aliasing effects by injecting periodic signals corresponding to any 
minima identified. 



4. Analysis of Individual Variable Sources 



4.1. J0939-2448 



J0939-2448 is an unusually blue T8 dwarf, exhibiting a K-suppressed s pectrum that 



likely arises from enhanced H2 a bsorption in a high-surfa ce gra vity atmosphere (jTinnev et al. 



2OO5I : iBurgasser et al.l l2006al ) . iBurgasser et al.l (l2008bl ) and iLeggett et al.l (120091 ) noted a 
mismatch between the near- infrared and mid- infrared spectrum of J0939-2448, and evidence 
of the source being overluminous, leading both studies to suspect that this source is an 
unresolved brown dwarf binary. Our observations of J0939-2448 over 2 hours on 2010 March 
28 (UT) reveal significant variability at K' with a semi-amplitude of 0.31 mags {S = 3.3), 
but no significant variation at J. The K' variability was the largest detected in our sample, 
but is tempered by the relatively low S/N 4) of the photometry over all dither positions 
for the entire exposure. Hence, we characterize this source as a marginal variable requiring 
confirmation. 



4.2. J1416+1348A 



J1416+1348A is an L6 dwarf flBowler et al.ll2010l : ISchmidt et al.ll2010r ) w ith a T7.5 spec- 
tral type binary companion at a projected sep aration of approximately 75 AU ( IBurgasser et al, 
2010bl : IScholzlboiol : iBurningham et al.lboioh . J1416+1348A is unusually blue in its NIR col- 
ors, possibly due to high surface gravity, low c ondensate opacities, unusually thin cloud s with 
large grain size, and/or subsolar metallicities (ICushing et al.ll2011t iDupuy fc Liull2012l ). Our 
observations of this source for 5 nights revealed signs of J-band variability only on the final 



night, with a semi-amphtude of 0.05 mags, a of 2.67 but with an S of 2.35, making this 
a marginal case. 



4.3. J1711+2232 



J1711+2232 is an optically-classified L6.5 dwarf ( Kirkpatrick et al. 200oh. for which 
near-infrared spectrosc opy has revealed signatures of CH4 absorption ( jNakajima et al. 2004; 
Burgasser et al.ll2010al ). This unusual absorption was interpreted in the latter study as an 
indication that this s ource is a spectral binary with L5 plus T5.5 components (it is unresolved 
in HST observations; iGizis et al.ll2003l ). Our short sequence of observations (0.7 hr) on 2010 
March 23 (UT) reveal signficant variability at J-band with a semi-amplitude of 0.10 mags [S 
= 2.5). The source also appears quite variable at K'; however, the comparison field sources 
are insufficiently stable to confirm this. Given the short timescale of the observation, we did 
not search for periodicity. 



4.4. J2139+0220 

J2139-I-0220 is a near-infrared classified T1.5 dwarf with unusually red near-infrared 
colors for its spectral type. The spectrum of J2139+0220 reve als significant CH^ abso rption 



at 1.4 fim, with diminished absorption in the H and K-bands iBurgasser et al.l (j2010al ). 



Our observations over three nights in October/November 2010 reveal persistent variabil- 
ity at J-band, with semi-amplitudes of 0.062 mag, 0.067 mag and 0.038 mag (5* = 3.6, 2.0 and 
2.2, respectively). Combining data from all three nights, we find a overall semi-amplitude of 
0.067 mag {S = 3.9). This was the most significant variation in our sample. 

Given its significant variation over a relatively long baseline, we performed both peri- 
odogram and PDM analyses on our multi-day light curve. Figure [9] displays the results of 
the former, for which we find significant peaks {la above significance of 1) at periods of 
11.7±0.5 hr, 7.42±0.22 hr, 5.44±0.17 hr and 4.26±0.12 hr. Figure [10] shows the light curve 
folded over 11.7 hr, driven largely by the overlap in data from October 28 and November 
1 toward the end of the phase cycle. Figure [TT] shows the light curve folded over 7.42 hr, 
driven by the overlap of one data point on November 1 with data on October 27. Figure 
[T2] shows the light curve folded over 5.44 hr, which has the most overlap in the night-to- 
night observations. Relative to the this last period, the others are found to be close to 2:1, 
4:3 and 3:4 harmonics. PDM analysis (Figure [T^ also shows similar significant periods at 
22.51±0.054 hr, 11.36±0.10 hr, 7.64±0.06 hr and 5.65±0.05 hr, which have ratios near 4:1, 
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2:1 and 4:3 about the shortest of these periods. 



The structure of the periodogram peaks reveals harmonics from peak-to-peak but also 
reveals double peaks for each significant point. Similarly in the PDM, the troughs are tripled 
for each significant point. In both cases, this is due to the sampling of the data, as the peaks 
themselves are not harmonics of each other. The separation between the periodogram peaks 
is approximately 0.3 hours for the shorter period peaks, to approximately 2.5 hours for the 
longer period peaks. The difference in the troughs of the PDM for each triple is between 
0.5 hours for the shorter periods, to 1 hour for the longer periods, which indicate sampling 
issues as observations were spaced out every 30 minutes. 

J2139-f-0220 has be en independently identified as a relatively large-amplitude variable 
by lRadigan et al.l ( 120121 ) based on several days of continuous near-infrared monitoring. That 
study finds a semi-amplitude of 26% at J-band — the largest detected to date — with a precise 
periodicity of 7.721±0.005 hours. This period is statistically consistent with the 7.42 hr 
and 7.64 hr periods found in our data, and Figure [H] compares our data with that from 
Radigan et al.l ( 12012| ) phased to 7.721 hr. We found a phase offset between our data and the 
Radigan et al. study of = 0.39. 

The Radigan et al. study also finds their light-curve to vary both in amplitude night- 
to-night and phase over several periods, which they interpret as arising from evolving cloud 
features. The amplitude of our variations of 0.067 mag is considerably smaller than the 0.15 
mag amplitude reported in the Radigan study. To assess if this difference is true variability 
or instrumentation and filter effects, we calculated the difference in relative brightness of 
this sources between the two filters used. We modeled the Gemini J-band filter as a box 



function with boundaries of 1.11 to 1.39 /im. According to the iRadigan et al.l (120121 ) study, 
the WIRC J-band filter was similar to the Mauna Kea Observatory (MKO) J-band filters, 

ved w ith the best-fit combination of cloudy and 



which we used here. These filters were convo 



clear atmospheres from ISaumon fc Marleyl (120081 ) that Radigan et. al found to reproduce 



the variability amplitudes at different wavelengths. The amplitudes were found to be the 
same for both filters indicating that the difference in the observed amplitudes is possibly due 
to intrinsic variability and changes in the cloud structure. 

Interpretation of these data must take into account the possibility that J2 139-1-0220 is 
a tight binary system since irradiation effects, magnetic interaction, or other orbital inter- 
actions could drive surface asymmetries. We tested this hypothesis with our NIRSPEC ob- 
servations. We cross-correlated each of the three NIRSPEC spectra wit h equivalent data for 



the T2 dwarf and radial velo city standard SDSSp J125453. 90-012247.4 (iLeggett et al.ll2000 



Zapatero Osorio et al.ll2007l ) following the methodology of iBurgasser et al. ( l2012h . The re- 
sulting velocities are listed in Table |5] and plotted in Figure [151 These data are marginally 
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consistent with a constant velocity of — 25.6±0.7 km s~^ {xf. = 1.83). A slight trend of 
— 5.4±2.0 km s~^ yr^^ also fits the data, but this is not statistically significant. Small varia- 
tions in the observed radial velocity of this source could arise from the same surface features 
ca using the light variatio ns. Assuming a single spot covering 26% of the surface as estimated 
by iRadigan et al.l (120121). radial velocity v ariations of up to 3 km s~^ could be discerned if 
V sin i fa 10 km s~^ ( ISaar fc Donahuelll997l ). In any case, we can definitively rule out the kind 
of short-period binary system that would give rise to interaction effects. Assuming that this 
interaction is magnetic in origin, so that the system is composed of two m = 0.05 MQbrown 
dwarfs se parated by at m ost twice their magnetosphere radii (a ~60 Rj^p = 4.3x10^ km ^ 
0.03 AU; ISchrijverl (j2009[ )) we can determine the maximum semi-amplitude of radial velocity 
via: 



irGm 
2P 



Gm 



2a 1 



sin[i] 



(6) 



where we assume the eccentricity e = and the inclination i = 90°; this yields , we 
obtain Ki ^ 27 km s~^. Assuming a 4a limit on the observed radial velocity variability 
amplitude of ~ 3 km s~^, this impies a maximum inclination imax = 6.2° if the variations 
are caused by binary interactions. The geometric probability of alignment for such a small 
separation P{i) = 1 — cosz = 0.0058 is sufficiently small to be deemed improbable. Indeed the 
3 km s~^ variability limit argues against a near-equal mass binary out to separation of 0.5 AU 
at 90% confidence limit. The period for this outer separation limit is 6 days; if we consider 
a binary whose orbit is synched to a period of 7.6 hours (~ 0.0042 AU or 8.83 Rjup), 
Ki 72. 'b km s~^, implying a maximum inclination of 2.4°, and we determine a probability 
of alignment P{i) = 1 — cosi = 0.00085. 

The apparently single nature of J2139-I-0220 does not necessarily contradict its identi- 
fication as a spectral binary candidate, however, as the dual spectral appearance may arise 
from distinct cloudy and non- cloudy regions on its ph otosphere. It is possible that other 
spectral binaries identified in iBurgasser et al.l (|2010aJ) may turn out to be similarly vari- 
able single L/T transition dwarfs with patchy atmospheres; 2MASS 1711+2232 may be an 
example of such a source. 



5. Discussion 

We have identified 4 variable L and T dwarfs in a sample of 15. This implies a rate of 
variability of 27^]^% with la errors based on binomial statistics. Tableland Figure [16] show 
how this variability rate compares to other "large" samples of 9 targets or more. Although 
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the studies were conducted with differen t filter s and report different sensitivity hmits, with 
the exception of iBailer- Jones fc Mundtl (120011 ) all agree with a variability fraction of 25- 
40%. For the 78 brown dwarfs in these 7 studies, 27 sources were variable, giving an overal l 
vari ability fraction of 35=l=5% . The >70% rates reported by iBailer- Jones &: Mundtl (1200 ll ) 
and iRockenfeller et al.l (120061 ) fro m a sample of 10 L dwarfs are considerably higher and in 
contradiction with the results of iGelino et al.l (120021 ). who observed a larger sample to the 



same depth and with a similar filter but found a variability rate of only 39~^iq%. In addition, 
we have examined the variability rate for distinct spectral ranges (Table E]), and find that 
the variability rate for each subset agrees with the overall rate of variability. We therefore 
conclude that the high variability fraction is likely spurious, and that existing data show no 
clear evidence of a change in the variability rate as a function of spectral type, even for the 
L/T transition. Nevertheless, it r emains true that th e largest amplitude variables detected 
so far— SIMP J013656.5+093347 JArtigau et allboooh and J2139+0220— are early-T dwarfs 
which show evidence of significant cloud structure. 

Improving on current statistics is challenging as current samples of known brown dwarfs 
are mostly located at high Galactic latitudes, where the number of field stars diminishes, 
making differential photometric analysis more difficult. As seen in Table |31 the number of 
comparison stars we were able to use was typically less than 5 and in some cases only 1 or 2. 
Coupled with the inability to detect small levels of variability (0.01 magnitudes and below), 
it is clear that infrared variability studies could be improved with larger field-of-view instru- 
ments. An instrument with an increase in linear size of 3 (> 9') would in principle increase the 
number of reference stars 10-fold. To our knowledge, there are no suc h instruments available 
on sm all (< 3m) tel escopes. Instrumen ts such as Flar ningos, UKIRT (iHurnphries fc Purkins 
1980h . NEWFIRM jAutrv et al. l l2003h . and IRIS 2 Jcillingham fc JonesI boooh on smaller 
telescopes would allow for more dedicated surveys of variable brown dwarfs. 



6. Summciry 

We have identified 4 variable L and T dwarfs in a sample of 15 observed: J1416-I-1348A 
, Jl 71 1+2232, and J2139+0220 in the J-band; and J0939-2448 in the K-band. Two of the 
variable sources were spectral binary candidates, and our detection of J2139+0220 confirms 
both the variability and period of this source as recently reported by Radigan et al. (2012), 
although we find that its variability amplitude has declined since its initial detection. We also 
find that this source has a stable radial velocity, suggesting that its spectral binary nature 
may derive from a patchy atmosphere rather than two distinct sources. We find a variability 
rate of 27l^Yo for our study and after examining 7 other similar studies (including our own). 
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we report an overall percentage of variability of 35 ± 5% for L and T type brown dwarfs. 
The challenge in observing low levels of variability necessitates the use of large field-of-view 
instruments on smaller and more dedicated telescopes. 
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Fig. 1.— NIRSPEC order 33 spectra for J2 139+0220 from (top to bottom) 2011 September 
7, 2011 July 6 and 2010 November 26 (UT). Spectra are corrected for barycentric motion, 
normalized and offset by integer constants. 



Fig. 2.— Master Gemini image of J2139+0220 for data taken 2010 October 27. The field of 
view is 2^9x2^9 with North up and East toward the left. J2139+0220 is indicated by the red 
circle, while comparison stars for differential photometry are marked in blue squares. The 
four ghost images around each source are due to the imperfect median combination of the 
overall sky frame. 
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Fig. 3. — Normalized differential magnitudes (differential magnitudes minus mean over an 
observing sequence) versus airmass (top) and seeing (bottom) for all observations and nights 
in J-band (left) and K-band (right). L-type dwarfs are plotted as diamonds while T-dwarfs 
are plotted as squares. We find no significant trends with these parameters for either spectral 
class. 
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Fig. 4. — Monte- Carlo simulation of variability for a signal-to-noise ratio S = 3. Two 
histograms representing a non- variable signal (unshaded) and a variable signal (gold shaded) 
are over plotted. The dashed line represents the top 1% of the non-variable signal and 
the dash dot lines represent the bottom 1% and 10% of a variable signal (left and right 
respectively) . 
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Fig. 5. — J-band light curves for all sources observed during the March 2010 run. Source 
name and observation date are listed to the right of each curve; uncertainties are overplotted 
on the data points in gray. For sources visited on multiple nights, light curves are displayed 
in the same color as the first night. Median differential magnitudes of the comparison stars 
for each source are indicated by the dotted lines 



-20- 



/\ I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



' 28 October 



'01 November 



J2148 + 40, 28 October 



J 1 November 



.J2224-01, 27 October 



I 28 October 



1 01 November 



I I I I 

' / 6 A [ 

Time (Hours) 



Fig. 6. — Same as figure [5] for tlie October 2010 run. 
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Fig. 7. — Same as figure [5] but ii'-band light curves for tlie Marcli 2010 run. 



- 22 - 




50 1D0 ^r^O 200 

Tins (mins) 



Fig. 8. — Light curve of J2139+0220. The blue hne at the top represents the total flux 
of the comparison stars, and is used to understand the general change in the atmospheric 
transparency throughout the course of the night. The uncertainties and the reduced chi- 
squared are listed. 
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Fig. 9. — J-band Lomb-scargle periodogram analysis of J2139+0220 spanning 3 days. The 
top panel shows the mean (red line) and uncertainty (green bars) of the raw periodogram 
based on the Monte Carlo simulation (all included periodograms are overplotted in light 
blue). The middle panel shows the window function for the 95% level. The bottom panel 
shows the significance of the mean periodogram, a significance of 1 is marked by the red 
dashed line. 
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Fig. 10. — Folded J-band light curve of J2139+0220 for all 3 days. The light curve from 
the 27th of October is shown in blue, the 28th of October shown in gold, and the 1st of 
November shown in pink. The top figure shows differential flux over time, while the bottom 
figure shows the differential magnitude over time. The light curve was folded over 11.7 hours. 
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Fig. 11. — Same as figure ITOl 
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Fig. 12. — Same as figure [TOl 
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Fig. 13. — Phase Dispersion Minimization results for J2139+0220. Periods in hours are 
plotted on the x-axis and ©, the ratio of variances in plotted on the y-axis. Minimum values 
represent true periods. 
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Fig. 14. — Folded flux light curve of J2139+0220 with data from Radigan et. al overplotted 
in blue with the + symbol, and our data in red circles phased to 7.72 hours, with a phase 
offset of 0.388, and a flux scaling ratio of 0.46 applied. 
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Fig. 15. — Radial velocity of J2139+0220 versus Modified Julian Date, spanning November 
2010 to September 2011. The reduced of the velocities, as well as the average radial 
velocity and radial velocity trend slope are listed. A constant -25.6 km s~^ radial velocity is 
the best fit case. 
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Fig. 16. — Variability rates over spectral type range for the studies listed in Table El We 
also indicate the mean and standard deviation for each spectral type range as stars. 
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Table 1. Observational Sample 



Source Coordinates Optical NIR 2MASS Type'' Refer. 

SpT SpT J Ks J-Ks 



J0141+1804 


01 


41 


03 


20 


+18:04:50.20 


LI 




L4.5 


13 


88+0 


03 


12 


49+0.03 


1.38+0.04 


Red 


1 


J0247-1631 


02 


47 


49 


90 


-16:31:12.60 






T2±1.5 


17 


19+0 


18 


15 


62+0.19 


1.57+0.27 


SpB 


2 


J0351+4810 


03 


51 


04 


37 


+48:10:46.80 






Tl±1.5 


16 


47+0 


13 


14 


99+0.12 


1.47+0.18 


Blue 


2 


J0407+1546 


04 


07 


07 


52 


+15:46:45.70 


L3 


5 




15 


48+0 


06 


13 


56+0.04 


1.92+0.07 


Red 


3 


J0805+4812 


08 


05 


31 


84 


+48:12:33.00 


L4 




L9+1.5 


14 


73+0 


04 


13 


44+0.04 


1.29+0.06 


SpB 


4 


J0825+2115 


08 


25 


19 


60 


+21:15:52.10 


L7 


5 


L6 


15 


10+0 


03 


13 


03+0.03 


2.07+0.04 


Red 


5 


J0939-2448 


09 


39 


35 


48 


-24:48:27.90 






T8 


15 


98+0 


11 




>16.6 


-0.62 ± 0.11 


SpB 


6 


J1106+2754 


11 


06 


11 


91 


+27:54:21.50 






T2.5 


14 


82+0 


04 


13 


80+0.05 


1.02+0.07 


SpB 


7 


J1315-2649 


13 


15 


30 


90 


-26:49:51.30 


L5 


5 




15 


20+0 


05 


13 


46+0.04 


1.73+0.07 


SpB 


8 


J1416+1348 


14 


16 


24 


09 


+13:48:26.70 


L6 




L6p±2 


13 


15+0 


03 


12 


11+0.02 


1.03+0.03 


Blue 


9 


J1511+0607 


15 


11 


14 


66 


+06:07:42.90 






T0±2 


16 


02+0 


08 


14 


54+0.10 


1.47+0.13 


SpB 


2 


J1711+2232 


17 


11 


45 


70 


+22:32:04.40 


L6 


5 




17 


09+0 


18 


14 


73+0.10 


2.36+0.20 


SpB 


5 


J2139+0220 


21 


39 


26 


76 


+02:20:22.60 






T1.5 


15 


26+0 


05 


13 


58+0.05 


1.68+0.07 


SpB 


3 


J2148+4003 


21 


48 


16 


28 


+40:03:59.30 


L6 




L6p 


14 


15+0 


03 


11 


77+0.02 


2.38+0.04 


Red 


10 


J2224-0158 


22 


24 


43 


80 


-01:58:52.10 


L4 


5 


L3.5 


14 


07+0 


03 


12 


02+0.02 


2.05+0.04 


Rod 


5 



''Red = unusually red; Blue = unusually blue; SpB = spectral binary. 



Ref eren ces. — flllWilson et al.l lIpOOSll f2')IChiu et al.l ll2006h fS 'l Eeid et alj||2008l') f4')lHawlev et al.l ||2002|') f5'llKi 
l l200(t) (6) iTinnev et all ||2005|) ('7') lLooper et al.l | |2007^ (8) iHalll ll2002h (9') lBowler et al.l 1 I2OIOI) (loT lLooper et al. 



•k £atrick et al.l 
ll2008bl) . 
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Table 2. List of sources and observation information 



Source 


Obs. 


No. of 


Time 


Avg Seeing 


Airmass 




Date (UT) 


dithers 


(hrs) 


(") 




J0805+4812 


24-Mar-2010 


25 


1.991 


1.61 


1.02-1.05 




28-Mar-2010 


22 


1.995 


1.44 


1.02-1.06 


J0825+2115 


23-Mar-2010 


3 


0.25 


4.04 


1.08-1.1 




27-Mar-2010 


6 


0.586 


1.97 


1.15-1.2 


J0939-2448 


28-Mar-2010 


20 


1.991 


1.42 


2.17-3.1 


J1106+2754 


23-Mar-2010 


18 


1.752 


3.66 


1.03-1.2 




29-Mar-2010 


21 


1.895 


1.69 


1.04 


J1315-2649 


24-Mar-2010 


6 


0.552 


1.73 


2.3-2.4 




27-Mar-2010 


6 


0.528 


1.96 


2.35-2.51 




28-Mar-2010 


8 


0.709 


1.89 


2.29-2.33 




zy-iviar-zuiu 





U.ooi 


Z.U4 




J1416+1348 


23-Mar-2010 


4 


0.369 


2.27 


1.11-1.13 




24-Mar-2010 


6 


0.562 


1.36 


1.13-1.17 




27-Mar-2010 


6 


0.52 


1.39 


1.1-1.12 




28-Mar-2010 


6 


0.518 


2.31 


1.12-1.15 




29-Mar-2010 


6 


0.526 


1.61 


1.19 


J1511+0607 


27-Mar-2010 


6 


0.523 


1.66 


1.29-1.37 


J1711+2232 


23-Mar-2010 


8 


0.683 


3.21 


1.03-1.04 


J2139+0220 


27-Oct-2010 


3 


1.603 


2.57 


1.25-1.9 




28-Oct-2010 


8 


2.671 


1.7 


1.24-1.72 




l-Nov-2010 


3 


0.216 


1.77 


1.22-1.25 


J2148+4003 


28-Oct-2010 


3 


1.227 


1.52 


1.01-1.12 




l-Nov-2010 


3 


1.248 


1.51 


1-1.2 


J2224-0158 


27-Oct-2010 


9 


2.565 


2.87 


1.29-1.77 




28-Oct-2010 


4 


2.449 


1.96 


1.3-1.52 




l-Nov-2010 


3 


1.217 


1.6 


1.29-1.37 
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Table 3. Results - J band 



Soiirce 


Obs. 
date. 


No. of 
data pts 


A 


var 


^med 


J-band 

s 


^comp 


X^ 




94 A/Tar 9010 




n 

u 


090 
uzu 





016 


1.25 


0.40 


Q 
O 





189 




98 A/Tar 9010 
Zo iVicli -ZU-LU 


22 


u 




Q 


022 


1.55 


0.78 


4 





474 




CombiiiGci 


A7 


n 


VJO-L 





015 


2.05 


1.08 


2 





737 




9^-Mar 901 


3 


Q 


015 





015 


1.00 


1.30 


3 





246 




97 A/T»r 901 


a 
o 


n 
u 


01 f\ 





018 


0.86 


0.41 


Q 
O 





470 




Combined 


Q 


u 







022 


1.51 


0.95 


Q 
O 





668 




98 A/T»r 9010 
Z o iVl dl Z U J- U 




n 
u 


141 





052 


2.71 


0.78 


Q 
O 





084 




9*^ A/Tar 901 
Z O- J-Vi dl - ZU i u 


1 s 


n 


022 





014 


1.57 


0.88 


Q 
O 





516 




9Q A/Tor 901 


91 


u 


099 





015 


1.47 


0.70 


Q 
O 


1 


678 




CombiriGcl 


oy 


n 







014 


2.87 


1.21 


Q 
O 





841 


J1315-2649 


24-Mar-2010 


6 





027 





021 


1.29 


1.00 


5 





448 




27-Mar-2010 


6 





Oil 





019 


0.58 


0.19 


7 





429 




28-Mar-2010 


8 





016 





020 


0.80 


0.45 


5 





549 




29-Mar-2010 


6 





095 





059 


1.62 


2.07 


3 





797 




Combined (24-28) 


20 





035 





022 


1.56 


0.43 


4 





898 


J1416+1348 


23-Mar-2010 


4 





006 





016 


0.38 


0.11 


4 





648 




24-Mar-2010 


6 





069 





114 


0.61 


0.24 


11 





378 




27-Mar-2010 


6 





033 





073 


0.45 


0.09 


10 





461 




28-Mar-2010 


6 





028 





025 


1.12 


0.78 


2 





738 




29-Mar-2010 


6 





054 





023 


2.35 


2.67 


2 





808 




Combined 


28 





052 





020 


2.54 


1.03 


4 


1 


097 


J1511+0607 


27-Mar-2010 


6 





036 





047 


0.77 


0.30 


10 





384 


J1711+2232 


23-Mar-2010 


7 





103 





041 


2.51 


1.81 


5 





387 


J2139H-0220 


27-Oct-2010 


8 





062 





017 


3.65 


8.72 


3 


1 


412 




28-Oct-2010 


3 





067 





034 


1.97 


3.80 


3 





883 




l-Nov-2010 


3 





038 





017 


2.24 


4.98 


3 





520 




Combined 


14 





067 





020 


3.35 


4.80 


3 





999 


J2148+4003 


28-Oct-2010 


3 





042 





038 


1.11 


0.94 


7 





382 




l-Nov-2010 


3 





010 





019 


0.53 


0.29 


14 





274 


J2224-0158 


27-Oct-2010 


8 





020 





Oil 


1.82 


1.10 


2 


2 


388 




28-Oct-2010 


3 





044 





020 


2.20 


6.04 


1 


5 


974 




l-Nov-2010 


2 





Oil 





015 


0.73 


1.12 


3 





234 
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Table 4. Results - K band 



Soijrce 


Obs. 


No. of 






K- 


band 








date. 


data pts 






S 




^comp 


comp 


J0805+4812 


24-Mar-2010 


25 


0.026 


0.019 


1.37 


0.60 


3 


0.80 




28-Mai-2010 


22 


0.083 


0.038 


2.18 


1.24 


4 


0.74 




Combined 


47 


0.068 


0.019 


3.52 


1.93 


2 


1.60 


J0825+2115 


23-Mar-2010 


3 


0.006 


0.015 


0.40 


0.18 


3 


0.06 




27-Mar-2010 


6 


0.020 


0.026 


0.77 


0.29 


4 


0.38 




Combined 


9 


0.012 


0.013 


0.90 


0.30 


3 


0.45 


J0939-2448 


28-Mar-2010 


20 


0.310 


0.094 


3.30 


2.51 


3 


0.73 


J1106+2754 


23-Mar-2010 


18 


0.031 


0.021 


1.48 


0.79 


3 


0.93 




29-Mar-2010 


21 


0.047 


0.021 


2.24 


1.18 


3 


1.06 




Combined 


39 


0.038 


0.02 


1.9 


0.75 


3 


0.96 


J1315-2649 


24-Mar-2010 


6 


0.025 


0.028 


0.89 


0.50 


7 


0.41 




27-Mar-2010 


6 


0.020 


0.021 


0.95 


0.70 


6 


0.40 




28-Mar-2010 


8 


0.100 


0.021 


4.76 


6.10 


2 


3.37 




29-Mar-2010 


6 


0.158 


0.028 


5.64 


12.1 


2 


8.71 




Combined (24-28) 


20 


0.078 


0.019 


4.07 


2.58 


3 


2.29 


J1416+1348 


24-Mar-2010 


6 


0.045 


0.036 


1.25 


0.66 


2 


0.85 




27-Mar-2010 


6 


0.022 


0.025 


0.88 


0.43 


2 


0.94 




28-Mar-2010 


6 


0.024 


0.046 


0.52 


0.25 


3 


1.24 




29-Mar-2010 


6 


0.045 


0.033 


1.36 


0.73 


2 


1.26 




Combined 


24 


0.052 


0.034 


1.55 


0.59 


2 


1.17 


J1511+0607 


27-Mar-2010 


6 


0.061 


0.036 


1.69 


1.84 


5 


0.47 


J1711+2232 


23-Mar-2010 


7 


0.593 


0.083 


7.14 


18.3 


4 


3.39 


J2139+0220 


27-Oct-2010 


8 


0.065 


0.034 


1.91 


1.91 


2 


3.65 




28-Oct-2010 


3 


0.035 


0.068 


0.51 


0.26 


2 


1.23 




l-Nov-2010 


3 


0.031 


0.044 


0.70 


0.89 


2 


1.01 




Combined 


14 


0.08 


0.04 


2.02 


1.21 


2 


2.09 


J2224-0158 


27-Oct-2010 


8 


0.063 


0.025 


2.52 


2.66 


1 


3.33 




28-Oct-2010 


3 


0.105 


0.022 


4.77 


15.8 


2 


5.09 




l-Nov-2010 


2 


0.012 


0.017 


0.71 


1.07 


3 


0.09 



Table 5. Radial Velocity Measurements for J2 139+0220 



Date 
(UT) 



MJD* 



RV 
(km s-i) 



2010 Nov 26 55526.197 -23.3±1.1 

2011 Jul 6 55748.543 -26.9±1.5 
2011 Sep 7 55811.377 -27.5±1.2 



''Modified Julian Date: Julian Date - 

2450000 



- 35 - 



Table 6. Variability fraction of L and T dwarfs for distinct spectral ranges 



Paper 


SoTvDe 
Range 


Total 
Observed 


Number 
Variable 


Fraction 


Detection Limit 
(mags) 


Filters 


Bailer- Jones and Mundt (2001) 


L0-L5 


10 


7 


70tl7% 


0.01 


I 


Gelino et al. ('2002') 


L0-L4.5 


18 


7 




0.01 


Ic 


Koen (2003) 


L0-T2.5 


12 


3 


25l^«% 


0.016 


Ic 


Enoch et al. f2003) 


L2-T5 


9 


3 


331^^% 


0.36 


Ks 


Koen et al. f2004~) 


L0-T8 


18 


6 


33tl^% 


0.02 


J,H,K 


Koen (20051 


L1-T6 


12 


3 


25l^«% 


0.02 


Ic 


Khandrika et. al (2012) 


L3.5-T8 


15 


4 


27i^'*% 


0.05 


J,K' 


Total 


L0-L5 


49 


18 








Total 


L6-T0 


13 


4 a 


301^-*% 






Total 


T1-T5 


11 ^ 


4 a 


36li5% 






Total 


T6-T8 


5 


1 


20tr^% 







^The total value is less than the sum as it contains no repeatedly-observed sources 
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